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Abstract. In this paper we describe a 1-dimensional variational ap- 
proach to the analytical construction of equivariant biharmonic maps. 
Our goal is to provide a direct method which enables analysts to com- 
pute directly the analytical conditions which guarantee biharmonicity 
in the presence of suitable symmetries. In the second part of our work, 
we illustrate and discuss some examples. In particular, we obtain a 1- 
dimensional stability result, and also show that biharmonic maps do not 
satisfy the classical maximum principle proved by Sampson for harmonic 
maps. 
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1. Introduction 

Harmonic maps are critical points of the energy functional 



where tp : (M, g) — > (N, h) is a smooth map between two Riemannian mani- 
folds M and N. In analytical terms, the condition of harmonicity is equivalent 
to the fact that the map ip is a solution of the Euler-Lagrange equation as- 
sociated to the energy functional i.e. 



The left member of (|1.2[) is a vector field along the map ip, or, equivalently, 
a section of the pull-back bundle <p~ x (TN): it is called tension field and 
denoted r (<p). In local charts, the tension field, which is the trace of the 
second fundamental form, is given by the following expression: 
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(1.1) 



trace \?dip = . 



(1.2) 



t j (v) = 9 t3 (V(dy»)) 
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(1.3) 
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where 

IVWJjy gxigxj l tJ + L a/3 Qx . . 

In (|1.4p . A/ r and N Y denote respectively the Christoffel symbols of the Levi- 
Civita connections of (M,g) and (N,h); note also that Einstein's convention 
of sum over repeated indices is adopted. In particular, inspection of (|1.4|) 
reveals that the harmonicity equation (|1.2p is a second order, elliptic system 
of partial differential equations. 

Every vector field v along <p defines a variation of ip by setting 

cp t (p) = exp^) (tv(p)) , t€ (-£,£) (1.5) 
(note that ipo — <£>)• If v is compactly supported, then 

dE(fp t ) 



V v E{ip) 



dt 



[r(ip), v)dv g . (1.6) 

t=0 JM 



In particular, it follows that p. 21) is equivalent to the vanishing of the direc- 
tional derivative in (|1.6j) for all w. Moreover, we point out that (|1.6j) implies 
the fact that the tension field r(<p) provides the direction in which the energy 
decreases more rapidly. 

The study of harmonic maps is a very wide area of research, involving a rich 
interplay of geometry, analysis and topology. We refer to [TJ |H1 EJ EH] for 
notation and background on harmonic maps and to [7] for a more recent 
bibliography. 

A related topic of growing interest deals with the study of the so-called bihar- 
monic maps: these maps, which provide a natural generalisation of harmonic 
maps, are the critical points of the bienergy functional (as suggested by Eells- 
Lemaire [10] ) 

Ei{f) = \ [ Wv)\ 2 dv g . (1.7) 

In [12] Jiang derived the first variation and the second variation formulas 
for the bienergy. In particular, he showed that the Euler-Lagrange equation 
associated to E 2 ((p) is 

r 2 (ip) = -J (t(<p)) = -Ar(ip) - traccR N (dv?, r((p))d<p = , (1.8) 

where J is (formally) the Jacobi operator of tp, A is the rough Laplacian 
defined on sections of tp^ 1 (TN) and 

R N (X, Y) = VxVy - VyVx - V[x,y] (1.9) 

is the curvature operator on (A, h). 

In this context, the biharmonic version of f| 1 .€>[) is (see [12]) 

dE 2 ( Vt ) 



V v E 2 (Lp) 



dt 



[T 2 (ip),v)dv g . (1.10) 
t=0 JM 



Therefore, ()1.8|) represents the vanishing of the directional derivative in (jl.lOp 
for all v. In particular, the bitension field T 2 (tp) provides the direction in which 
the bienergy decreases more rapidly. 
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We point out that (|1.8[) is a fourth order semi-linear elliptic system of differen- 
tial equations. We also note that any harmonic map is an absolute minimum 
of the bienergy, and so it is trivially biharmonic; thus, a general working plan 
is to study the existence of biharmonic maps which are not harmonic that we 
shall call proper biharmonic. We refer to [13] for existence results and general 
properties of biharmonic maps. 

Fourth order differential equations are of great importance in various fields. 
By way of example, we cite an instance which is well-known to civil engineers: 
the structural problem of a beam resting on an elastic soil is amenable to 
the following differential equation, in which the unknown function Y(x,i) 
represents the response of the beam in position x, at the time t: 

Q2 / Q2 Y \ Qly 

3^{ EI[X) ^) +m ^ + kY = fM ' (LU) 
where EI(x) measures the flexural rigidity (Young modulus), m and k are 
two constants depending respectively on the material of the beam and on the 
elasticity of the soil; the right member f(x,t) represents the external load 
(see [1] for details). 

In general, we encounter enormous difficulties to study fourth order differ- 
ential equations. In particular, the presently known instances of biharmonic 
maps have always been obtained essentially by means of geometric intuition 
and simplification (see, for example, [HH]). 

In this order of ideas, equivariant theory deals with special families of maps 
having enough symmetries to guarantee that harmonicity reduces to the study 
of a second order ordinary differential equation (we refer to [8j I15j for back- 
ground and examples). 

In this paper we shall explain how this general framework can be adapted to 
include the study of biharmonic maps. Although the study of the resulting 
fourth order ordinary differential equation remains, in general, a very difficult 
task, one of our aims is to provide analysts with a direct approach which per- 
mits them the computation of the relevant equations avoiding, in particular, 
the often heavy burden of dealing with an advanced ricmannian geometric 
machinery. As an application, we obtain a 1-dimcnsional stability result and 
also show that proper biharmonic maps in general do not satisfy Sampson's 
maximum principle. 

Acknowledgement. The authors wish to thank the referee for comments and 
suggestions that have improved the paper. 

2. The 1-dimensional variational approach 

In order to introduce the variational context in which we work, let us first 
recall some basic facts concerning equivariant maps and their associated 1- 
dimcnsional variational problem. Various, rather general approaches are pos- 
sible (see, for example, jU [5J [T5] ) : here we adopt the setting of [15], which is 
largely sufficient for our purposes and, at the same time, will ease our task 
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to make this article as self-contained as possible. In particular, we shall con- 
sider maps / which are equivariant with respect to Riemannian submersions: 
keeping, just for this short introduction to equivariance, the notation of [T5] , 
that amounts to require that the following be a commutative diagram: 

Ei — - — y E'2 

TrJ (2.1) 

Mi — ^— >• M 2 . 

where 7Tj : Ei — > Mi , i = 1,2 are Riemannian submersions. We say that 
a vector field v along / is basic if: for any p £ E\, v(p) is horizontal with 
respect to 7T2 and there exists a vector field v along / such that 

dn 2 {v{p))=v{-K X (p)) . (2.2) 

We need the following adaptation of a result of Xin (see Theorem 6.5 in [15] ) 
to the biharmonic case: 

Lemma 2.1. Let f : E\ — > Ei be an equivariant map with respect to Rie- 
mannian submersions. If its bitension field t 2 is a basic vector field, then f 
is biharmonic if and only if it is a critical point of the bienergy with respect 
to all equivariant variations. 

Proof. Although the proof follows [T5] step by step, we report it here for 
the sake of completeness. We only have to prove that if / is a critical point 
with respect to all equivariant variations then / is biharmonic. Let h be any 
equivariant function on E 2 with compact support and h its induced function 
on M 2 . Consider the following variational maps 

ft(p) =exp f(p) (th(f(p))T 2 (f)) , t£ (-£,£). 

By construction the corresponding variational vector field is (h o f)r 2 {f). 
We now show that /* is an equivariant variation of /. Since T 2 (f) is basic, 
there exists a vector field v(f) along / such that d-K 2 {T 2 {f){p)) = v{f){-K\{p)). 
Thus, we have 

7T2°/ f (p) = CXp n2of{p) (dlT 2 (th(f(p))T 2 (f))) 

= cxp />l(p)) (t/io /(7Tl(p))C(/)(7ri(p))) 

= MMp)), 

with 

f t (q) = cxp f - {q) (tho f(q)v(f)(q)) , q £ M± , t € (-£, e) , 

where h is such that h(p) = h(jr 2 (y)), for any y £ E 2 . Now, by using the first 
variational formula (|1.10j) . we obtain 

dE 2 (f t ) 



= 



dt 



Ei 



(r 2 (/), (hof) T2 (f))dv g 
(hof)\r 2 (f)\ 2 dv g , 
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which implies, by the assumptions, that T%(f) = 0. □ 

In order to describe our applications, we restrict our attention to the case that 
the base manifolds in (|2.ip are 1-dimensional. More specifically, we consider 
equivariant maps ip a : (M, g) — > (N, h) and denote by 

a:[a,b]^K. (2.3) 

the associated map between the base manifolds. Note that this function a 
may have to satisfy suitable boundary conditions dictated by the geometry of 
the problem. If these maps have enough symmetries (in particular, to ensure 
that the tension field is basic) (see [5] for details, and $3] for examples), then 
the energy functional (jl.ll) takes the following form: 

E v (a) = ( L(t,a,a)dt , (2.4) 

J a 

where a denotes derivative with respect to t, and L(-, •, •) is a differentiable 
function depending on the geometry of the problem under consideration. 
Now, according to (the harmonic maps version of) Lemma l2.1[ the harmoni- 
city of the map ip a is equivalent to the fact that a is a critical point of the 
so-called reduced energy functional (|2.4p . That is to say, the function a must 
be a solution of the Euler-Lagrange equation associated with (|2.4p , which is 
the following: 

dL d f ° L ^ . (2.5) 



da dt \ da / 

For future reference, it is useful to recall how equation (|2.5[) is derived from 
(|2.4I) . The first step in this direction is to recognize that the requirement that 
a be critical is equivalent to the vanishing of the directional derivative 



— [Er(a + hp)]\ h=0 

for all compactly supported (differentiable) variations /3 : 
we compute (|2.6p explicitly: 



d 

dh 



[E*(a + hp)]\ h=0 = 




h/3)dt 



hp)\ h= 



(2.6) 
Next, 



h=0 



dt 



where, in order to obtain the fourth equality of (|2.7[) 
that, since (3 is compactly supported: 



(2.7) 



we have used the fact 
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By way of summary, we conclude from (|2.7p that the vanishing of the direc- 
tional derivative (j2.6p for all compactly supported variations /? is equivalent 
to the Euler-Lagrange equation (|2.5D . 



We are now in the right position to extend this setting to the framework of 
biharmonic maps. In the equivariant context, the bienergy functional E2(f) 
introduced in (| 1 . T[) takes the following form: 

E$(a) = I L(t, a, a, a) dt , (2.9) 

J a 

where, as above, the function a : [a,b] —¥ M. may have to satisfy suitable 
boundary conditions, and L(-, •, •, •) is a differentiable function depending on 
the geometry of the problem under consideration. Now, according to Lemma 
12. 1[ the condition of biharmonicity for an equivariant map (p a with basic 
bitension field is equivalent to a being a critical point of the reduced bienergy 
functional (12.911. Next, we can state our most useful result in this context: 



Theorem 2.2. A differentiable function a : [a, b] — > R is a critical point of the 
reduced bienergy functional (|2.9j) if and only if it is a solution of the following 
differential equation: 

A / PiT \ ^2 / P.T \ 

= . (2.10) 

Proof. In (|2.8|) we have already noticed that, if f3 is a compactly supported 
function, then 

rb dL • f b ( d dL" 



dL 


d 




d 2 




da 


~dt 


Kda) 


+ dt 2 


(da) 



da^ dt --] a {dtda-)P dt - ^ 
In a similar vein, we shall also need the following equality, which holds when 
both f3 and f3 are compactly supported: 

b dL ■■ f b f d 2 dL 



m^ dt -JA^da) Pdt - (2 ' 12) 

Now, we proceed to the verification of (|2.12|) : we start from 



from which we deduce: 



b dL ■■ , f b ( d dL 



d^ dt =-L [dtda-^ dt - ^ 



Next, we observe that 







b d_ Wd_ dL\ fl ] M _ f b fd 2 dL\ n „ 
dt. \di~da~) 



ft 



dt =l ^mF dt + J a (dtm^ dt - (2 - 15) 



By using (|2.15p in (|2.14p we easily obtain (|2.12p . We are now in the position 
to prove the theorem: a function a is a critical point of the functional E 1 ^ (a) 
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if and only if 



dh 



{E%(a + hp)]\ h=0 = 



(2.16) 



for all (diffcrcntiablc) variations (3 : [a, b] — > R, with (3 and $ compactly 
supported. Next, we compute the directional derivatives in (|2.16|) explicitly: 



dh 



[E%(a + hp)]\ h=0 = 




L(t, a + h/3,a + h/3,a + h(3) dt 



h=0 



— L(t,a + hP,a + hp,a + hP)\h=a) dt 



8L 
da 
'dL 
da 



0- 



dL • dL n 

P + TTTT P) dt 



da da. 
d dL d 2 dL 
dt da dt 2 da 



dt , (2.17) 



where, in order to obtain the fourth equality of (|2.17p . we have used both 
(|2~TT|) and ([2~T2"|) . Now it follows that the validity of (|2"JB)) for all (3 is equi- 
valent to the fact a is a solution of (|2.10p . as required. 

□ 

By way of conclusion, an equivariant map ip a with basic bitension field is 
biharmonic if and only if a is a solution of (f2~T0|) . 

Remark 2.3. Our approach require explicitly that the equivariant map ip a 
have basic bitension field. Therefore, it would be interesting to obtain a com- 
plete geometric characterization of the situations in which this property is 
satisfied. In general, this appears to be a rather technical and difficult pro- 
blem. However, here we point out that there are important, large families 
of equivariant maps for which it is immediate to conclude that tension and 
bitension field are both basic. In particular, that occurs when the Rieman- 
nian submersions Tti : Ei — > Mi , i = 1, 2 in (|2.1j) are determined by isometric 
actions of Lie groups, and the equivariant maps <p a , when restricted to the 
fibres endowed with the induced metric, are harmonic maps with constant 
energy density. All the examples in below are of this type. 

Remark 2.4. For biharmonic curves in a Riemannian manifold, the Euler- 
Lagrange method was investigated in [5] where, with a different method, an 
equation of type (|2.10|) was derived. 

Remark 2.5. In some important geometric applications, such as the study of 
biharmonic immersions, for instance, it is necessary to consider a variant of 
(|2.9j) . in which the unknown function a is replaced by a curve. More precisely, 
we have to consider: 

,6 

E% ( a i) = / L(t,ctj,dj,ctj) dt , j = l...r. (2-18) 
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In this case, the argument given in the proof of Theorem 12.21 applies again 
and leads us to the conclusion that the critical points of the functional (|2.18|) 
are precisely the solutions of the following system of ordinary differential 
equations: 



dL 


d 


(—) 


d 2 


(—) 


da 3 


~ dt 


{da J 


+ dt 2 


\dajj 



3. Equi variant biharmonic maps and applications 

In this section we discuss some examples and their applications to problems 
concerning stability and maximum principle. 

Example 3.1. Equivariant maps from the flat 2-torus T 2 to the 2-sphere S 2 . 

We write the flat 2-torus as a product manifold: 

T 2 = (S 1 x S 1 ,dj 2 + d6 2 ) , 0<7,6><2tt. (3.1) 

Next, we describe the 2-sphere S 2 by means of polar coordinates: 

S 2 = (S 1 x [0,tt], sin 2 a df3 2 + da 2 ) , 0</3<27r,0<a<7r. (3.2) 

We consider equivariant maps ip a : T 2 — > S 2 of the following form: 

( 7 , 0)H>(fc 7 , a(9)) , (3.3) 

where k S Z is a fixed integer, and a is a differentiable, periodic function to 
be determined (in this case, the period must be equal to 2tt). We also note an 
alternative way to describe an equivariant map of type (|3.3[) . More precisely, 
considering the canonical isometric embedding of S 2 into R 3 , we can rewrite 
(|3.3p as follows: 

(7, 0) H. ((sina(0)) (sinfc7), (sina(6>)) (cosfc 7 ), (cosa(8))) . (3.4) 

If one wishes to allow the unknown function a to take values outside the 
interval [0,7r], then the form (|3.4I) is to be preferred. In any case, a direct 
computation shows that the tension field of these equivariant maps can be 
written as follows: 

d 

T(w a ) = \a(0) - k 2 sina(0) cosa(0)l — . (3.5) 
L J da 

Therefore, writing a instead of a(8), we find that the reduced bienergy (|2.9|l 
associated to this family of equivariant maps is given, up to a constant, by: 

E$(a) = / [a-k 2 sina cosa] 2 dO (3.6) 
Jo 



k 4 

a 2 + — sin 2 (2a) - k 2 d sin (2a) 



d6 
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Next, we can compute the condition of biharmonicity by applying directly 
(12.101) to (|3.6p . We obtain the following fourth order ordinary differential 
equation for a: 



a (4) -d [2k 2 cos(2a)] + a 2 [2k 2 sin(2a)] + — sin(2a) cos(2a) = . (3.7) 

At present, as we have already pointed out in the introduction, we do not 
dispose of general methods to carry out a satisfactory qualitative study of 
this type of equations. However, in this example we can observe that there 
are some trivial solutions, namely: 



(i) a =£ - , where I = 0, 1, 2 



.... 7T 3 7T , . 

ii) a = — or a = — . 3.8 

v ' 4 4 



The solutions in (|3.8p (i) correspond to harmonic maps. By contrast, those in 
(|3.8|) fii) represent proper biharmonic maps. Although these biharmonic maps 
were known by other methods (indeed, they can be obtained by projection 
of T 2 onto S , followed by a biharmonic immersion (see |13j)). our approach 
leads us to a rather surprising 1-dimensional stability property. More pre- 
cisely, we can prove the following result: 

Theorem 3.2. Let ip a : T 2 — >• S 2 be a proper biharmonic map as in (|3.8j) (ii). 
Then a is a local minimum for the reduced bienergy functional (|3.6I) . 



Proof. We have to compute the second variation of the reduced bienergy 
functional (|3.6p at the point a = (tt/4) (which will be denoted by a*). We 
find: 

' {E%(a* + hV)]\ h=0 



V'E%(a*)(V,V) 



dh 2 



d^_ 

dh 2 



{hVf 



k -2/^ 

4 v 2 



2hV) 



dP_ 
dh 2 



-hVk 2 sin(- + 2/iV) 
h 2 V 2 + ^ cos 2 (2hV) 



dO 



h=0 



2n 



-hVk 2 cos (2ft y) 



/i=0 



2 V 2 



2 V 2 k 4 



de 



(3.9) 



It follows from the computations in (|3.9j) that the solution a* is strictly stable, 
as required to end the proof. The case a = (3n/4) is analogous, so we omit 
it. 

□ 



Remark 3.3. We point out that the reduced bienergy of the biharmonic maps 
of Theorem 13.21 is strictly positive. Therefore, they are local, but clearly not 
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absolute minima. This fact makes it reasonable to conjecture that Mountain 
Pass techniques may be used to prove the existence of (unstable) equivariant 
biharmonic maps of rank 2. 

Now we are going to consider a situation which represents a generalization 
of our previous Example 13.11 More precisely, we shall work with warped 
products of the following type: 

(S m x(a,b), f 2 (r)g sm +dr 2 ) , (3.10) 

where (S m , gs m ) denotes the Euclidean unit m-sphere, and / is a smooth, 
strictly positive function on the interval (a,b). For instance, if f(r) = r, 
then our warped product is (a subset of) the Euclidean space R m+1 . The 
cases f(r) — sinr and f(r) = sinhr correspond respectively to S m+1 and to 
hyperbolic space H m+1 , while f(r) = c > is a cylinder. 
We can introduce the following family of equivariant maps between warped 
products of type (|3.10|1 : 

Va : (S m x (a, b), f( r )g sm +dr 2 ) ->■ (S n x (c,d), h 2 (a) g S n + da 2 ) 

( 7 , r) ^ (* A ( 7 ), a(r)) , (3.11) 

where ^ is a so-called eigenmap of eigenvalue A . That means that ^ \ : 
S m — > S n is a harmonic map with constant energy density equal to (A/2). 
Important examples of eigenmaps are: the identity map of S m (A = m), the 
k-fold rotation e l9 e lke of S 1 (A = k 2 ); and, also, the Hopf fibrations 
S 3 -> S 2 , S 7 -> 5 4 and S* 15 -> 5 8 , with A equal to 8, 16 and 32 respectively. 
Now, the unknown function a(r) in ()3. 1 If) has to be determined in such a 
way that ip a be a biharmonic map. To this purpose, a calculation based on 
(|1.3p shows that the tension field of an equivariant map of the form p. lip is 
given by the following expression: 



•■ / \ f(r) . , v , h(a) h(a) 



Therefore, in this case the reduced bienergy (|2.9p is given, up to an irrelevant 
constant factor, by the following expression: 

f{r) ... x h(a)h(a)" 2 



J a 



a(r) + m d(r) — A 



/(r) v y / 2 (r) 



f m {r)dr. (3.13) 



Now, a direct application of (|2.10j) to p.!3p yields the biharmonicity equation 
for a. But, since this equation is in general rather long and complicated, we 
shall write it down and analyze in a few particular instances only. We start 
with: 

Example 3.4. Equivariant maps between Euclidean spaces (compare with [1]). 

This situation corresponds to the case that f(r) = r and h(a) = a in (|3.13l) . 
If we also perform the following change of variable: 

r = e\ teR, /3(t) = a(e t ), (3.14) 
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then the reduced bienergy functional (|3.13[) takes the following form: 

E$(P) = I fj9 + (m-l)4 -\fi\ 2 e (™- 3 )* dt . (3.15) 



By way of example, we assume that "J a : S 3 — > S 2 is the Hopf map, so 
that m = 3 and A = 8. A simple computation shows that the biharmonicity 
equation (|2.10p becomes: 

/3 (4) - 20/3 + 64/3 = , (3.16) 

which admits the following explicit representation of solutions: 

p(t) = ci e 4t + c 2 e- 4t + c 3 e 2t + c 4 e" 2t . (3.17) 

We observe that, when either c\ or C4 is different from zero, the map which 
is associated to (|3.17l) is proper biharmonic: this fact is actually part of a 
general principle (Almansi's property), which states that multiplication of a 
harmonic function on R m+1 by the factor r 2 generates a proper biharmonic 
function (see [13]). 

Example 3.5. Equivariant maps from cylinders to Euclidean spaces (violation 
of the maximum principle). 

Now we illustrate the case of equivariant maps ip a : S m xl-} K n+1 . This 
instance corresponds to the choice f(r) = 1, r € K and h(a) = a in (|3.13l) . 
so that the reduced bienergy functional (|3.13[) becomes: 

E%(a) = I [a-\a)] 2 dr. (3.18) 

Next, applying (|2.10|) to (|3.18() . we find that the biharmonicity equation for 
this case is: 

a (4) -2Aa + A 2 a = (3.19) 
(note the similarity between equation (|3.19l) and a time independent version 

of srm \ 

Equation (|3.19p admits the following explicit representation of solutions: 

a (r) = ci e^ r + c 2 e~^ r + c 3 r e^ r + c 4 r e~ VXr . (3.20) 

We observe that, when either C3 or C4 is different from zero, the map asso- 
ciated to (|3.20p is a proper biharmonic map of full rank (n + 1) (of course, 
provided that we are using an eigenmap of rank n) . 

In |14j , Sampson proved the following maximum principle for harmonic maps 
(to state this result, we keep the notation of [H]): 

Theorem 3.6. Assume that f : M — > Y is a harmonic map, with q = f(p). 
Let S be a C 2 hypersurface in Y passing through q, at which point we assume 
that the second fundamental form is definite. If f is not a constant mapping, 
then no neighbourhood of p is mapped entirely on the concave side of S . 
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In order to construct our counterexample, we consider an equivariant bihar- 
monic map arising from the following solution of p.20p (A > , ^\ of full 
rank): 

a(r) = (VAr) sinh(VAr) + eS^ r) , rel. (3.21) 

It is easy to verify that the solution a(r) in (|3.2ip admits a strictly positive 
absolute minimum point, say r — ro- Therefore, the image of our biharmonic 
map is entirely contained in the concave side of S = dB r0 (O) . This fact 
shows that proper biharmonic maps do not satisfy the Sampson maximum 
principle of Theorem 13.61 
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